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c§ : ABSTRACT 
^ ■ 

>■ ■ We analyze the local field of stellar tangential velocities for a sample of 42 339 

^ ■ non-binary Hipparcos stars with accurate parallaxes, using a vector spherical 

CN ' harmonic formalism. We derive simple relations between the parameters of 

the classical linear model (Ogorodnikov-Milne) of the local systemic field and 
Q ■ low-degree terms of the general vector harmonic decomposition. Taking ad- 

^ ' vantage of these relationships we determine the solar velocity with respect to 

■ the local stars of (Vx,VV,14) = (10.5, 18.5, 7.3) ± 0.1 kms"\ not corrected 

■ for the asymmetric drift with respect to the Local Standard of Rest (LSR). If 
j_| ' only stars more distant than 100 pc are considered, the peculiar solar motion 

■ — ' is {Vx,Vy,Vz) = (9.9, 15.6, 6.9) ± 0.2 kms~^. The adverse effects of harmonic 

leakage, which occurs between the refiex solar motion represented by the three 
electric vector harmonics in the velocity space and higher-degree harmonics in 
the proper motion space, are eliminated in our analysis by direct subtraction of 
the reflex solar velocity in its tangential components for each star. The Oort's 
parameters determined by a straightforward least-squares adjustment in vector 
spherical harmonics, are A = 14.0 ± 1.4, B = —13.1 ± 1.2, K = lA ± 1.8, and 
C = -2.9±1.4 kms"^kpc~\ The physical meaning and the implications of these 
parameters are discussed in the framework of a general linear model of the ve- 
locity field. We find a few statistically significant higher degree harmonic terms, 
which do not correspond to any parameters in the classical linear model. One of 
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them, a third-degree electric harmonic, is tentatively explained as the response 
to a negative linear gradient of rotation velocity with distance from the Galactic 
plane, which we estimate at ~ —20 kms~^kpc"^. A similar vertical gradient of 
rotation velocity has been detected for more distant stars representing the thick 
disk {z > 1 kpc), but here we surmise its existence in the thin disk at 2; < 200 
pc. The most unexpected and unexplained term within the Ogorodnikov-Milne 
model is the first-degree magnetic harmonic representing a rigid rotation of the 
stellar field about the axis —Y pointing opposite to the direction of rotation. This 
harmonic comes out with a statistically robust coefficient 6.2 ±0.9 kms^^kpc^^, 
and is also present in the velocity field of more distant stars. The ensuing up- 
ward vertical motion of stars in the general direction of the Galactic center and 
the downward motion in the anticenter direction are opposite to the vector field 
expected from the stationary Galactic warp model. 

Subject headings: astrometry — Galaxy: kinematics and dynamics 



1. Introduction 

The growing amount and quality of available astrometric data and the widening horizons 
of the astrometrically known Galaxy begin to call for more general and accurate models of 
the Galactic velocity field than the classical Oort and Lindblad approximation. The classical 
Oort's constants A and B representing the local angular velocity of rotation and its radial 
gradient have been determined many times on different data sets with increasing accuracy, 
but only the advent of the Hipparcos catalog made it possible to put this estimation on a 
systematically rigid, global reference frame related to a non-rotating extragalactic frame. 
A somewhat heuristic generalization of the Oort's formalism, known as the Ogorodnikov- 
Milne model (OMM hereafter), describes the local field of stellar velocities as a combination 
of all six possible shears and three rigid rotations, incorporating the A and B parameters, 
and adding seven more degrees of freedom to the fitting model. Not all of these additional 
parameters have a clear physical meaning in terms of Galactic dynamics. To our knowledge, 
only two additional coefficients [K and C) have been paid close attention to in the recent 
literature. The remaining five OMM parameters have been tacitly considered non-essential 
and presumably small. Furthermore, as a linear approximation to a complex and, on a large 
scale, nonlinear Galactic velocity field, the Oort's formalism and the OMM are meaningful 
only locally, within reasonably small distances from the Sun. As we expand the limits of 
precision astrometric catalogs, we i nevitably encounter signifi cant deviations from the linear 



model. Following the suggestion of lOlling fc Merrifieldl (119981 ). one could consider the OMM 
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parameters to be functions of heliocentric distance, or, in a more general approach, of galactic 
cylindrical coordinates. Alternatively, as we investigate in this paper, one can use a more 
versatile formalism for describing the observed velocity field in order to put this problem on 
a more rigorous and systematic footing. 



The paper by IVityazev fc Shuksto I (120041 ) marked an important advance in the search 



for a better way to describe the systematic field of tangential velocities derived from Hip- 
parcos proper motions and parallaxes for tens of thousands stars. They employed vector 
spherical harmonic functions, which constitute an orthogonal basis of the space of contin- 
uous vector functions on unit sphere. When sampled over a large and sufficiently uniform 
set of points on the celestial sphere, the discretized vector harmonic functions remain nearly 
orthogonal, which makes the fitting algorithm the most stable and accurate in implementa- 
tion. Since the Oort's and OMM parameters turn out to be simply the fitting coefficients 
of the corresponding low-degree vector harmonic functions, the physical meaning and inter- 
pretation of the estimated quantities remain straightforward and simple. The mathematical 
formulation allows us not only to accurately estimate the uncertainties associated with the 
estimated parameters (arising mostly from the stochastic component of the velocity field), 
but also, for the first time, to investigate the covariances. But the main advantage of this 
method is that it is easily expandable to handle more complex and nonlinear fields. The aims 
of this paper are: 1) using the vector harmonic method, to re-determine the fundamental 
parameters of the local velocity field in a rigorous fashion, viz., the peculiar solar velocity and 
the differential rotation coefficients A and B; 2) to estimate other OMM parameters; 3) to 
find out if the observed velocity field bears statistically significant, higher-degree non-OMM 
harmonic terms that may reveal interesting dynamical phenomena, such as vertical rotation 
gradients. Galactic warp and large streams of stars. 



2. The peculiar velocity of the Sun 

The dynamic LSR is defined as a fictitious point currently at the position of the Sun 
in the Galactic plane, which moves along a perfectly circular orbit in a hypothetical ax- 
isymmetric potential. This definition is a mere theoretical concept, because the Galactic 
potential is not exactly axisymmetric, and there are no perfectly circular orbits. An alterna- 
tive empirical definition of the LSR is the average motion of a sufficiently large, sufficiently 
homogeneous and sufficiently dynamically mixed sample of stars centered on the Sun. These 
two definitions are quite different, and in fact, are contradictory with regard to the shape 
of the reference orbit. Indeed, the empirical average reference orbit is markedly eccentric 
due to the stellar age-dependent asymmetric drift (see Fig. [1]). The relative velocity of the 
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Sun with respect to the empirical LSR is directly derived from our heliocentric astrometric 
observations, whereas determination of t he solar peculiar velocity with respect to the hypo- 



thetical circular orbit is more involved (IMihalas fc Binneyl Il98ll : iDehnen fc Binney I Il998l ) 



and relies on additional astrophysical or dynamical considerations. 

If the Sun moves with a velocity Vq = (Vx, Vy, Vz) relatively to the average motion of 
the local stars, the heliocentric velocity field of these stars contains a streaming motion in the 
opposite direction, that is, —Vq. As far as tangential velocity components are concerned, this 
streaming motion is a dipolar vector field on the celestial sphere, and its exact representation 
via vector spherical harmonics (see Eqs. IA6I) is 

Vr = VxEl + VyE^' - VzEl (1) 

Thus, the signature of the solar motion is confined to the first three electric harmonics 
only. We will see in the subsequent paragraphs that the other fundamental parameters of 
the velocity field are represented by magnetic harmonics and electric harmonics of higher 
degree, so that this effect is clearly separated in the vector harmonic space. It is worth noting 
that Vq is estimated from tangential velocities in physical units, i.e., from yu/II, where ^ is 
the proper motion magnitude, and 11 is the parallax. A small admixture of halo stars and 
runaway stars with very high spatial velocities can perturb this determination. 

We start with selecting 42 487 stars from the main Hipparcos catalog with statistically 
robust parallaxes (n/cr(n) > 5) and without any indicators of binarity. To avoid extra-high 
velocity perturbers, we reject 148 stars with tangential velocities greater than 150 kms~^ in 
either galactic component. A set of 24 vector harmonic functions is then fitted to the global 
vector field of 42 339 stars by a direct least-squares solution. The solar velocity components 
are simply read from the fitted coefficients of the first three electric harmonics. The results 
of this estimation are specified in Table [H along with a sample of more distant stars with 
n > 10 mas (24 327 stars). 

The estimated velocity components Vx and Vz are fairly similar for the two samples, 
in dicating a negligible d epend ence on distance. They are also very close to the determination 



by iDehnen fc Binney I (Il998l ). The estimates of Vy (in the direction of Galactic rotation) 
are very different between the two samples, beyond the possibility of a statistical fluke. The 
more distant stars move faster with respect to the Sun than the stars closer in. In either case, 
the stars move slower i n this direction by n a ore th an 10 kms~^ than the circular motion of 



the LSR determined by IDehnen &: Binney I (119981 ). This very prominent effect is attributed 



to the asymmetric drift of nearby stars, discussed in more detail in the next paragraph. 
Generally, there are three different reasons of physical and technical kind for our estimation 
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of Vy to be biased: 



1. the asymmetric drift; 



2. the vertical gradient of rotational velocity T{z); 



3. the mixing of non-orthogonal harmonics. 



3. The physical meaning of Oort's constants 



As long as a relatively small local area of the Galaxy around the Sun is considered 
(r << Rq, where Rq is the galactocentric distance of the Sun), it is appropriate to expand the 
systemic velocity field of stars in a Taylor series over the galactocentric cyhndrical coordinates 
{p,6,z). Assuming that the local field is planar, that is, all systemic motions are in the 
galactic plane, the velocity vector v{p, 9, z) — a{p, 9)+u{p, 9, z), where the former component 
a is radial with respect to the Galactic center, and the latter, uj is tangential and orthogonal 
to the former. Note that the radial component is assumed to be independent of z, that is, 
that there is no radial systemic motion depending on the distance from the plane. Retaining 
only first-degree terms in the corresponding Taylor expansion, one can write 



where the subscript denotes the corresponding parameters at the Sun's location 
{po,9o,Zo). We left the dependence of the rotational velocity a; on z in its generic form, 
since on physical grounds, this dependence is expected to be symmetric around the plane, 
and can not be represented by a simple linear term. 

Retaining only terms to O(^), these model relations can be rewritten more conveniently 
in the heliocentric coordinates (r, i, b) introduced in Appendix A: 



a{p,9) 
u;(p,9,z) 



ao + S{p- po) + e sm{9 - 9q) 

uJo + a{p- po)+ P sm{9 - 9o) + T{z) 



(2) 



a{r) — ao — Sr cos £ cosb + e — sin ^ cos 6 

Po 



(3) 



u(r] = ujq — ar cosi cosh + P — sin £ cosh + T(z) 

Po 



r 



Apart from the reflex peculiar motion of the Sun treated in Section [2|, we observe the 
heliocentric velocity field Av = a + uj — Sq — i^o- Projections of this vector field onto the local 
tangential coordinate frames (r^, fh) introduced in Appendix A, are 



Av ■ Te = rcos2£cos6 r cosb + {uj{r) — ujq) cos£ (4) 

2po 2po 



CL\T] 

r sm2£ cos 6+ (a(r) — Oq) sin^ 



2pi 







Av ■ Tfc = r sin 2i cos b sin b — (cuif) — luq) sin i sin b 

2po 

r cos b sin b H r cos 2£ cos b sin 6 + {a{f) — oq) cos £ sin b 



2po 2po 

Substituting the model ([3]) into Eqs. H] and retaining only terms to 0{^), we obtain 
after some toil the general expansion 
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rH^ - T{z) El 



-1 



(5) 



2po 6 2po 

where we made use of the functional forms of the vector spherical harmonics in galactic 
coordinates specified in Appendix A. Disregarding for now the T{z) terms, let us compare 
this equation with the classical e xpansion of the velocity field via the fundamental Oort's 
constants (e.g.. lTorra et al.ll200Cll ). which in the vector harmonics notation takes the form 



(6) 
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The slope of the local rotational velocity curve is readily derived as a = —{A-\-B). 
Since most of the recent and Hipparcos-based estimations arrive at A ^ —B, the rotation 
curve is locally almost (but not exactly) flat. It is also usually adopted that the local angular 
rotation velocity 60 is just the difference of the constants A and B. In fact, however, 

eo = ^ = A-B + ^. (11) 
Po Po 

So, the local azimuthal shear of the radial motion e contributes to the constants A and 
B and affect the determination of the angular velocity of the Galaxy. Presumably, this shear 
is small (e << cjq), and our estimations are not hampered too much. The interpretation 
of the constants C and K is more complicated. Traditionally, the C constant is called the 
shear, and the K constant the local expansion (or dilation). In fact, we find that 



C + K = ^ (12) 
Po 

K-C = 6 + —. (13) 
Po 

Presumably, there is no outward or inward bulk motion of stars around the Sun (oq = 0), 
and the LSR, empirically defined as the mean motion of a large homogenous heliocentric sam- 
ple of stars, moves on a circular orbit. In numerous and somewhat conflicting determinations, 
it appears that, generally, C 7^ —K, which indicates a nonzero systemic eccentricity of the 
local field. It is especially important in this case to accurately estimate the uncertainties of 
estimation, which may be larger than the estimates, as discussed in Section |H The dilation 
is better characterized by the difference — C in Eq. [T31 where S can be interpreted as the 
radial heliocentric expansion, and jS/pQ as the azimuthal expansion. These two components 
can not be separated from proper motions alone. 



4. Parameters of the velocity field 

Using the vector harmonic formalism to describe the tangential velocity field on the 
celestial sphere (Appendix A) and the general expression for Ogorodnikov-Milne model 
(Eq. IBSp . makes the estimation of model parameters quite straightforward. The Hippar- 
cos proper motions for our initial set of 42 339 non-binary stars with accurate parallaxes are 
converted to tangential velocities in the galactic coordinate system. Each velocity vector 
generates two condition equations, one for the longitudinal component and the other for 
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the latitudinal component ti. The coefficients of the expansion IA2I are the unknowns of 
the condition equations, which are solved by the least-squares method. The main source 
of the solution uncertainty is the physical dispersion of individual velocity vectors related 
to peculiar orbital motions, since the astrometric errors of proper motions are small in the 
Hipparcos catalog, and the number of stars is large. The velocity dispersion is known to 
increase with age; it is also larger for thick disk stars than for thin disk stars. Instead of 
dealing with triaxial dispersion ellipsoids for various stellar populations, we take an empirical 
and robust approach to error estimation in this least-squares adjustment. We fit 24 vector 
harmonic functions up to degree 4 to the general vector field and consider this expansion 
to represent the systemic part of the velocity field. The residuals of the tangential velocity 
vectors represent the stochastic part of the field. Dispersions of velocities are computed 
from these residuals, separately for the fe and fb components, as half-differences between the 
0.84 and 0.16 quantiles on each distribution. These quantities substitute standard deviation 
parameters for the markedly non-Gaussian velocity distributions. The resulting dispersions 
are a{vi) = 24.0 kms~^, a{vb) = 16.1 kms~^. The condition equations in longitude and in 
latitude are weighted with these quantities, respectively. 

There are a few important technical notes to be made on this estimation problem. 
The solar peculiar velocity vector is determined directly from the tangential velocity field, 
expressed in units of kms~^ O, in which case only the ffist three electric harmonics 
are of essence. The Oort or Ogorodnikov-Milne parameters describe a velocity field which 
grows linearly with distance from the Sun, and has distance r in its functional form flB2p . 
The corresponding decomposition is done in the proper motion field, or, as we do it in this 
paper, the tangential velocities can be used in the observational part of the equations, but 
the harmonic functions are pre-multiplied with distances for each star. In the latter case, 
the distribution of sample stars on distance is taken into account automatically, and the 
harmonic coefficients have the desired dimension of kms~^ kpc~^. But before performing 
this distance- weighted least-squares estimation, the relative velocity of the Sun with respect 
to the stellar centroid should be subtracted for each star. This step proves to be of crucial 
importance because of the adverse effects of the harmonic leakage, discussed in § 14. 1[ 

The results of vector harmonic estimation are specified in Table |5] for the original set 
of 42 339 stars, and for 24 327 stars more distant than 100 pc. All harmonic coefficients 
corresponding to Ogorodnikov-Milne parameters are shown, as well as other statistically 
significant coefficients which have no counterparts in the linear model. By statistically 
significant we conservatively mean a quantity larger than its formal error multiplied by 2.5. 
Having the significance so defined we state that only three or four model OMM parameters 
are significant, and three extra non- linear parameters. The estimated parameters agree very 
well between the two sets, indicating no considerable dependence of the velocity field on 
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distance within this fairly small volume. The slope of the rotation curve, using the results 
for the larger sample, is a = —{A + B) = —1.0 ± 1.8 kms~^ kpc^^. This implies that the 
speed of rotation declines very slightly with galactocentric distance, but the conclusion is 
not reliable statistically. From Eq. [TT], ignoring the possible contamination by the azimuthal 
shear e, the local angular rota tion is A — B = 27.1 ± 1.8 kms~^ kpc~^, in fine agr eement 
with iFeast fc Whitelockl (119971 ). Assuming a distance po = 7.9 kpc for the Sun (see IVallee 
20051 . and references therein), the speed of rotation is ujq = 214 ± 14 kms~^. 



The K constant is insignificant for both samples, but the C constant is marginally 
significant, especially for more distant stars. The systemic outward motion oq = Po{C + K) 
appears to be negligible for the general sample, but the more distant stars seem to exhibit 
an inward motion of an = —42 = b 20 kms^^. This result is qualitatively consistent with 
the estimation by iHanson I (119871 ) who found a progressively smaller solar velocity toward 
the Galactic center with respect to stars at higher latitudes. If this inward motion in the 
outer part of the astrometrically known Galaxy is real, it may be somewhat counterbalanced 
by the small but persistent dilation (expansion) of the local stellar aggregate, at -ft" — C ^ 
4 ± 2 kms~ ^ kpc~^. Association s of young stars expand by virtue of their initial velocity 
dispersions ( jMakarov et al. 1120041 ). and the presence of the young Local Association could 
be invoked to explain the local expansion. It is worth emphasizing that the accuracy of the 
available astrometric data on the local stellar field is still insufficient to establish these subtle 
effects with certainty. In fact, the barely noticeable K and C constants may be related to 



the intermediate-scale streams of stars pe rmeating the so . 



the general pattern of Galactic rotation. iFamaey et al. I (120051 ) present a scrutiny of such 



ar ne ighborhood, rather than to 



streams or superclusters, based on the best available radial velocity and astrometry data for 
K and M giants, including the Hyades, Sirius, Hercules and B streams. Apart from strong 
evidence for asymmetric drift for evolved stars, they find, interestingly, that the centroid 
velocity of the Sun Vx with respect to giants is ~ 10 kms~^, in agreement with our present 
and other previous estimations, but it drops to only ~ 3 kms~^ when all the major streams 
are excluded. This difference r nay be interpre t ed as a net outward radial motion of the 



streams (see also their Table 1). iFamaey et al. I (120051 ) point out that the members of these 



streams have a spread of ages and other physical characteristics, and the streams must be 
dynamically induced. The authors raise the question of how the standard solar motion can 
actually be defined if the motion of even the oldest and supposedly dynamically mixed stellar 
populations is subject to unknown dynamical agents perturbing their orbits? We think that 
the stellar streams are legitimate parts of the local velocity field, and that it makes sense to 
define the velocity centroid and the solar peculiar motion in much the same way as it has 
been done before, keeping in mind that dynamical mixing and relaxation may be a mere 
theoretical idealization, as well as a circularly moving LSR. 
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4.1. Harmonic leakage 

As specified in Table [2l we find only three significant terms in the general vector har- 
monic decomposition beyond the Ogorodnikov-Milne model, viz., H2^, E^^, and Ef. All 
other estimated harmonics, including all third and fourth degree terms, are well below 2.5a. 
Thus, we find little evidence of nonlinear patterns in the motion of local stars. The actual 
velocity field progressively deviates from the linear approximation of the model with he- 
liocentric dista nce. Furthermore, the rota tion curve may have local wiggles and curvature. 



as discussed in lOlling fc Merrifieldl (119981 ). One way of tackling this problem is to build a 
more complex model in which the Oort's constants are actually functions of coordinates, to 
be determined from observations. We take a different approach in this paper, determining 
empirically a vector harmonic decomposition and trying to interpret those terms that appear 
to be statistically significant. 

Before embarking on analysis of the emerging nonlinear harmonics (and the unexpected 
linear term H^^), we should examine a technical, but crucial problem in the determination 
of model parameters. The stellar velocity field bears a strong signal in the classical terms 
representing the reflex solar motion and the Oort's constants A and B. These terms are rep- 
resented in our model by specific vector harmonics (Appendix B). The strong signal in the 
physically meaningful low-degree harmonics in the velocity space can leak into higher-degree 
vector harmonics in the proper motion space, resulting in spurious detections of nonlinear 
effects. This inevitably happens because the sampled vector harmonics are not independent 
for any inhomogeneous discrete set of points. This problem has two somewhat different 
aspects. In classical applications, when only proper motions are known from observation, 
the mean parallax of nearby stars varies across the sky because of the real clumps in num- 
ber density (the Gould Belt, large associations, spiral a rms), as well as th e non- uniform 



interstellar extinction. This difficulty was first spott ed by lEdmondson I (119371 ). and later in 



vestigated in more detail by lOUing fc Dehneru (120031 ). In the latter paper, a nice example is 



presented, how the longitudinal variation of the mean parallax, described by a Fourier series, 
makes the simple dipolar pattern of the solar motion to contribute to the terms that would 
be empirically defined as the A and B constants. Our analysis is free of this complication, 
because we use accurate trigonometric parallaxes from the Hipparcos catalog, and perform 
the estimations of the solar motion in the velocity space, and of model parameters separately 
in the proper motion space. But there is another, more basic reason to be concerned about 
the harmonic leakage. The lack of uniformity in the number density of stars on the sky itself 
makes the vector harmonics mutually dependent within either coordinate component. 

Mathematically, the problem can be viewed as a lack of orthogonality of the sampled 
harmonics. The degree of non-orthogonality is quantified by the correlation coefficients 
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readily computed from the off-diagonal elements of the covariance matrix. For our nearby 
stars, the largest physical effect is the solar motion expressed by the first-degree electric 
harmonics, and the cross-talk of these terms with other harmonics of higher degree may 
generate false positive detections. 

We set up a dedicated numerical experiment to prove that this contamination may hap- 
pen unless appropriate precautions are taken. We use the same general set of Hipparcos 
stars as before, but the actual observed proper motions are replaced with simulated vectors, 
computed from the reflex solar motion only, estimated in §2J A harmonic decomposition of 
the simulated velocity field produces the same velocity dipole in the first electric harmonics, 
and zero for the rest of harmonics, which only shows that the software works correctly. But 
when a similar decomposition is carried out in the space of distance-weighted harmonics, 
as described in §U a number of spurious terms emerge, viz., (significance level 8.80"), 
Hi (5.4a), (3.2a), E^'^ (8.9a), E^^ (4.0a), E^ (7.6a), E| (8.7a), and E| (3.8a). The 
appearance of the i?2^^ and E^^ harmonics is especially worrisome, because they may carry 
some physical information, as discussed in subsequent paragraphs. The simplest way to 
get rid of most of the harmonic leakage effect is to subtract the reflex solar motion from 
all tangential velocities prior to a model parameter fitting. Ideally this eliminates the per- 
turbations from the dominating dipolar terms. Note that existing correlations between the 
sampled distance- weighted vector harmonics that we use to determine OMM parameters, 
do not affect the results in a systematic way, because the least-squares solution is unbiased. 
The major adverse effect of these correlations is an enhanced propagation of random and 
possibly systematic errors from our observational data. 



5. Vertical gradient of rotational velocity 



The apparent relative velocity of the Sun in the direction of galactic rotation (Vyq) 
varies with distanc e of reference f ield stars from the Galactic plane. This remarkable fact 
was establish ed by Hanson I ( 19891) from proper motions of fairly distant stars, and recently 
confirmed by Girard et al. ( 20061 ). who used absolute proper motions of giant stars in the 
direction of south Galactic pole. The thick disk dominates between z = 1 and 3 kpc, where 
the rotational lag of field stars is found to follow a nearly linear dependence on vertical 
height, accompanied, predictably, by a growth of velocity dispersion in the X direction. The 
slope of the lag, from both cited papers, is estimated at —30 kms~^kpc~^. Girard et al. 
also offer a dynamical interpretation of this phenomenon, finding it consistent with a general 
model of the Galactic potential. The sample of nearby Hipparcos stars considered in this 
paper is practically limited to 200 pc, and is dominated by thin disk stars. Is there a similar 
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vertical gradient of rotational velocity for the thin disk? 

Evidently from Eq. IA2[ a vertical lag affects the determination of the centroid velocity 
Vy expressed by the dipole vector harmonic Ei^, because T{z) is negative everywhere except 
z = 0. If the velocity of rotation falls off with increasing z, the relative solar velocity Vyq 
should grow with distance due to the admixture of high-z stars. This is not what we find 
in Table [H where the more distant stars (11 < 10 mas) appear to rotate faster than the 
overall sample of stars. However, a more accurate consideration reveals that for a number 
of possible functional forms of T{z) (e.g., F ■ \z\, T ■ l-zp), the most characteristic response is 
expected in the E^^ harmonic, because the E^^ harmonic is too sensitive to the choice of 
centroid solar motion. Our choice of velocity in Table [H consistent with the estimation for 
the more distant half of Hipparcos stars, is justified by the fact that the OMM parameters 
are determined in the distance-weighted (or proper motion) space where distant stars are 
more significant, and whatever kinematics anomalies the nearest stars may have, has little 
bearing on the OMM estimation problem. The rotation gradient dipole E^^ emerges with 
a robust positive coefficient e^^ = 11.25 ± 1.17 kms~^kpc~^, which is consistent with a 
negative gradient of T{z). 

We performed direct simulations of the vector harmonic response to a linear gradient 
r(z) = r ■ |2;| for different values of F and the height of the Sun above the plane Zq. A 
good match with observations was found for F = —20 kms~^ kpc~^and Zq = 15 pc, which 
yielded a set of coefficients e^^ = 11.41 ± 1.17, = -2.04 ± 0.47, and ^ = 0.75 ± 0.22 
kms~^kpc~^, all the rest 21 harmonics being insignificant. Both electric harmonics are in 
good agreement with our fit for all star, whereas the coefficient is fairly close to the fit 
(—1.09 ± 0.47, not shown in Table [2]). Therefore, a linear gradient of rotational velocity of 
the thin disk of roughly —20 kms~^ kpc~^is a plausible explanation to the corresponding set 
of vector harmonic terms beyond the Ogorodnikov-Milne model. The detected pattern of 
tangential velocities of Hipparcos stars consistent with this interpretation is shown in Fig. [2J 

6. Warp and the origin of H^^ and 

The Milky Way disk is warped, as has been established from the distribution of stars and 
neutral hydrogen. In this respect, our Galaxy is not different from many other spiral galaxies 
exhibiting a range of warp distortions. The origin of galactic warps is not clear; a number 
of hypotheses have been proposed, including the tidal interaction of the disk with the dark 
matter halo, the influence of the bar, and the perturbation from a major satellite galaxy. 
The Sun appears to be close to the line of nodes of the Milky Way warp, and the upper 
rim of the disk is at ^ = 90° (in the rotation direction). The height of the midsection above 
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the plane is quadratic with galactocentric distance in the model of iDrimnael et al.l (l2000f). 
w(p) = (p — pwY/15 kpc for p > pyj, and zero for p < p^. According to iMomany et al 



(l2006l ). the warp begins well within the solar circle (p^ < P0), and the line of nodes deviates 



from the solar radius by 15°. 

The single most unexpected result of our analysis is the strong model parameter L13 
(Table [2]), represented by the coefficient of the first-degree magnetic harmonic -fff^, that is, 
a rigid rotation around the direction —Y (see Eqs. IA7p . The sign of this parameter implies 
that the stars move upward in the direction of the Galactic center, downward in the opposite 
direction, away from the center at the north pole, and toward the center at the south pole. 
The signal-to- noise ratio on this parameter is about 6. The extra statistically significant 
term h^^ detected by us in the local velocity field may be related to the former. The pattern 
of tangential velocities generated by these two magnetic harmonics, Q.21rH^^ — 1.20 ri?2^^, 
is shown in Fig. [31 The main effect of the higher degree harmonic is that the axis of rotation 
lies below the plane at roughly h = —20°, nearly obliterating the motion in the north pole 
region, but retaining the galactocentric motion near the south pole. The most conspicuous 
features are the general upward motion of stars in the direction of the galactic center, and 
the downward motion of stars at £ = 180°. It is tempting to relate these two unexpected 
components to a kinematic signature of the Galactic warp. 

The shape of the warp, as traced by the distribution of neutral hydrogen and dust, 
implies that the stars in the solar region are involved in a general upward motion, since the 
starting rim of the warp is within the solar circle (p^ = 6.5 kpc). This common motion is 
indistinguishable from the vertical solar refiex motion, but the model also implies a radial 
gradient of the upward velocity, Vz,warp (p — p^), which is detectable in the proper motion 
field. In the near-plane zone, the differential warp motion manifests itself as a downward 
stream at £ ~ 0, and an upward stream at i ^ 180°. Obviously, the pattern in Fig. [3] is 
completely inconsistent with this prediction. Assuming for simplicity that the Sun lies on the 
line of nodes, the linear Taylor expansion of the local velocity field ([2]) should be expanded 
to include a vertical [Z) component of velocity. 



A(p) = Ao + A(p - po) ^ Ao - A r cos i cos b. (14) 
The corresponding tangential velocity field is 

Atv = -Ar cosicos'^bTb = -Ar{^H{^ - ^El). (15) 

Thus, the differential warp motion is expressed by two OMM parameters, L13 and 
M13. These two fitted parameters yield discrepant estimates of the warp velocity gradient. 
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A = —2h^^ ~ —12 kms~^kpc~^, and A = 663 — +4 kms~^kpc~^. The former estimate 
from the magnetic harmonic has the wrong sign, and its modulus is too large for a credible 
differential warp. The electric harmonic has the right sign, but it nearly vanishes for more 
distant stars (11 < 10 mas). Thus, the kinematical model of Galactic warp does not furnish 
an adequate explanation to the presence of magnetic harmonics H^^ and i?2^^. Samples of 
larger volumes are needed to find out if these two harmonics are not a local feature, a nd to 
find evidence of warp in the motion of field stars. Interestingly, iDrimmel et al.l (120001 ) also 
found a negative vertical motion of distant OB stars in the direction of Galactic anticenter, 
in obvious contradiction to the predicted warp motion. A non-stationary warp is one of the 
possibilities considered by them. A precessing line of nodes is conceivable, but we find it 
difficult to reconcile the observed pattern of vertical motion, should it bear on the subject at 
all, with a plausible precession model. It appears instead, that the line of nodes is stationary, 
but the shape of warp changes to its opposite every 50 Myr or so, curling this way and the 
other. 



7. Conclusions and Back to Astrometry 

Hipparcos stars with accurate trigonometric parallaxes represent only a tiny fraction of 
the Galactic population. Half of stars considered in this paper are within 112 pc, and 75 % 
are within 160 pc. The narrow horizon of our selection limits the accuracy of vector harmonic 
terms describing the local tangential velocity field in the most general and systematic fashion. 
Only several major kinematical parameters can be determined with confidence from such a 
limited data set. We determine the relative solar velocity with respect to all stars in our 
selection, and to stars with measured distances greater than 100 pc. The latter determination 
yields Vq = (9.9, 15.6, 6.9) ± 0.2 kms~^, which we add to the tangential velocities of all field 
stars before performing a general decomposition of the velocity field onto vector spherical 
harmonics. This decomposition provides values of A = 13.8 ± 1.4 and B = —13.4 ± 1.2 
kms~^kpc~^ for the fundamental Oort's constants of differential galactic rotation. Since 
A + B ^ 0, the local rotation velocity curve is nearly flat. Assuming a galactocentric distance 
of 7.5 kpc for the Sun, a rotation velocity of luq = 204 ■ [po/7.5 kpc] kms^^ is derived. 

Among other linear OMM parameters, we detect, most unexpectedly, a strong signal 
carried by the first-degree magnetic harmonic H{^, which describes a rigid rotation of the 
stellar field around the axis —Y opposite to the direction of galactic rotation. The estimated 
rate of this rotation is roughly 6 kms~^ kpc~^, or 1.3 masyr"^ in proper motions along the 
principal galactic meridian. Another unexpected magnetic harmonic, i?2^^, nearly cancels 
out the outward motion at the North pole predicated by the former harmonic, but retains 
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the strong inward motion around the South pole, and the counter vertical motions near 
the Galactic plane. Differential vertical velocities naturally arise from a kinematic model 
of the Galactic warp, but we find that the sign of the local rotation is opposite to what is 
required to raise the rim of the Galaxy above the plane in the first and second quadrants. In 
other words, the local stars are expected to move upwards due to the warp, but we detect a 
negative differential rotation. Analysis of velocity fields in a much larger volume of space is 
needed to make sure that this discordant rotation is not a local feature, which would have 
crucial consequences for our understanding of the physics of the warp. 

Only three statistically significant vector harmonic terms beyond the Ogorodnikov-Milne 
model are detected in this paper. One of them, the electric multipole ^3"^, is of special note, 
since, together with a positive residual dipole E^^ in the direction of galactic rotation, it is 
likely to advertise a vertical gradient of rotation velocity. A similar gradient of rotational lag 
of —30 kms^^kpc^^ has been found in proper motions of more distant stars representing 
the thick disk population, but never reported for the thin disk dominating our sample. We 
estimate a gradient of ^ —20 kms~^ kpc~^ for our sample of nearby stars limited to 200 pc. 
This result requires verification on a larger sample of thin disk stars extending to 1 kpc. 

Our concluding remark is that estimation of subtle effects in the local kinematics per- 
taining to the Galactic structure and formation history is based on the assumption that the 
Hipparcos proper motion data is free of large-scale systematic errors at > 1 masyr"^. None 
such errors have been reported in the literature, which is not a strong argument because 
Hipparcos remains unparalleled at its level of globa l astrometric accuracy. The raajor c ata- 
logs of proper motions Tycho-2 (jUrban et al. Il2000l ) and UCAC (jZacharias et al. II2004I ) are 
calibrated on Hipparcos stars; therefore, systematic distortions of Hipparcos astrometry, if 
any, are just copied over to these catalogs. Radio astrometric observations with VLBI have 
recently advanced to a comparable level of accuracy in positions and proper motions, and 
being di rectly tied to the IC RF, provide an independent test for the Hipparcos reference 
system (IBoboltz et al. II2006I ). This important external check is unfortunately limited by 
the small number of optically bright radio stars, but the available accuracy of VLBI proper 
motions (approximately 1.7 masyr"^) enables Boboltz et al. to state that the relative spin 
of the Hipparcos proper motion system is much less than 1 masyr"^ about each axis. This 
result confirms that the strong magnetic harmonic H^^ representing a spin around the —Y 
direction, is not an artefact. Another significant astrometric development of late is the SPM3 
catalog, which provides high quality absolute proper m otions for a large sa mple of distant 
and faint stars, albeit in a small fraction of the sky (IGirard et al. II2004I ). This catalog 
provides an independent view of the local stellar velocity field in the surveyed area of the 
sky. 
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A. Vector spherical harmonic decomposition of a proper motion field 

As customary in studies of Galactic dynamics, we make use of the Galactic coordinate 
system (X, F, Z) in which the X axis is pointing toward the Galactic center, the Y axis 
toward the direction of Galactic rotation, and the Z axis toward the north pole. For each 
star, a triad of unit vectors (f, 77, 7^) is defined, with 



( cos ^ cos 6 \ / —sm£ \ / — cos £ sin 6 \ 

sin £ cos 6 = cos£ ^ \ — sin ^ sin 6 , (Al) 

sin6 / \ / \ cos6 / 

where and define the tangential coordinate directions toward increasing galactic 
longitude i and the north pole, respectively. The proper motion vector of the object is 
traditionally projected onto the locally tangential coordinate vectors, that is, jl = cos 67^+ 

A global proper motion field of a large set of celestial objects can be represented by the 
expansion 



00 n 



(A2) 



n=l m=—n 



where £ and b are galactic longitudes and latitudes, and E!^ are orthogonal vector 
harmonics which we call magnetic and electric vector harmonics respectively. These vector 
harmonics are derived via partial derivatives of the scalar spherical harmonics over angular 
coordinates, viz.: 



H^{i,b) 



dS^{£,b)^ 1 dS^{£,b) 



db 



-n 



cos b dl 



1 asrM)J 

" H — n 



cos b d£ 



db 



Spherical harmonics S^^ are counted by degrees n — 0, 1, . . . and orders m — 
1, . . . ,n. Explicitly, 



(A3) 



-n, —n + 



2n + l{n-m)\ 
S"^ = W (-P™(sm6) cosm£, m > 0, 



2tt {n-\-m)\ 
P'^nb) 



2n + 1 „o 



47r 



m = 0, 



(A4) 
(A5) 



where are the associated Legendre polynomials. The first pair of vector harmonics are 
generated by from the scalar zonal harmonic , with the electric component cos brb and the 
magnetic component cosfer^. The electric vector harmonics for n <= 2 in angular coordinates 
are 





= — cos i fi + sin ism bfb 




El 


= cos b ff, 




El 


= smife + cos £ sin 6 7^ 






= 6 cos 2i cos bfe — 6 sin 2i cos b sin b 


n 




= — 3 cos i sin bfi — 3 sin £ cos 2b fh 




4° 


= 3 cos b sin bfh 




El 


= 3 sin i sin bfi — 3 cos i cos 2b fl 




El 


= —6 sin 2i cos bfi — 6 cos 2i cos b sin 


b 



and the magnetic vector harmonics for n <= 2 are 

= sin £ sin 6 + cos £ t% 

Hi = cos b fi 

HI = cos £ sin bfi — sin £ f^ 

= —6 sin 2£ cos b sin bfe — 6 cos 2£ cos b 

= —3 sin £ cos 26 7^ + 3 cos £ sin 6 fi 

i?2 = 3 cos b sin 6 7^ 

Hi = —3 cos £ cos 26 f^ — 3 sin £ sin b fb 

Hi = — 6 cos2£ cos6 sin6ff + 6 sin 2£ cos^T^, (A7) 



Note that the normalization coefficients of the generic spherical harmonics flA4p are 
omitted in these formulas. 
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B. Ogorodnikov-Milne model via vector harmonics formalism 



The linear Og orodnikov-Milne model of the local velocity field in its general matrix form 
can be written as (Idu Mont 1119771 ) 



Mu 


Mi2 


Mi3 ^, 




f 


—L12 


-L 




M22 


M23 









-L 


Mi3 


M23 


M33 ) 




K Ll3 


L23 






V = { M12 M22 M23 f + L^2 -L23 f, (Bl) 




where V is the systemic part of the velocity field as a function of the position vector r. 
It is assumed in the following that positions are determined with respect to the solar sys- 
tem barycenter; any shift of the coordinate system origin results in an additional constant 
translation term. For convenience, the matrix of transformation is split into the symmetric 
(shear) part M and the antisymmetric traceless (rotation) part L. It is readily seen that 
the matrix M describes the gradient-type distortions of the field, and the L part represents 
rigid rotations, or spins, around the three coordinate axes. 

After a small manipulation, the tangential velocity components are 

Vi/r = {V ■ Tt,)/r = (-Mil + M22)sin£cos£cos6 + Mi2Cos2£cos6 

— Mi3 sin i sin b + M23 cos i sin b + L12 cos b + L13 sin i sin b — L23 cos i sin b 
Vb/r = {V ■ Tb)/r = (—Mil cos £- M22 sin^£ + M33) sin b cos b — M12 sin 2i sin b cos b 

+M13 cos i cos 26 + M23 sin i cos 26 + L13 cos i + L23 sin i 

Comparing these equations with the trigonometric expressions for low-degree vector spherical 
harmonics (Appendix A), the following relations of proportionality are established 

M12 oc e^^ 

M23 (X -62^ 

Mi3 oc -el 

Di oc el 

D2 oc el 

L12 oc h\ 

Li3 oc h^^ 

L23 oc -hi, (B2) 

where ^ ^ 

D, = -(Mil - M22); D2 = M33 + -(Mil + M22). (B3) 
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Clearly, all nine parameters of the Ogorodnikov-Milne model can not be determined 
from a proper motion field, since an isotropic dilation (Mn = M22 = M33) results in radial 
velocities only. This is why only eight pa rameters of the model appe ar in the vector harmonic 
decomposition of a proper motion field (IVityazev fc Shuksto 1 120041 ) . 

We further establish the relations between the constants of the Ogorodnikov-Milne 
model and the four constants {A, B, C, K) of the Port' s two-dimensional model by matching 
the terms in the proper motion equation of the latter (ITorra et al.ll2000l ) 



4.741 /i = {Acos2i cost + B cosb - C sm2£ cost) Ti 

+ {—Asin2l£smbcosb — C cos 2£ sin 6 cos 6 — sin 6 cos 6) Tf,. (B4) 

The corresponding equation for the more general Ogorodnikov-Milne model in terms of vector 
harmonics is 

4.741 /i = Li^H^' + Li2H'i ~ L2sHl 

+ ^M^2E2^ - ^M2^E^^ - ^D2E^2 - Im.sEI + ^D,El (B5) 

where Mu = A, L12 = B, Di = C and D2 = K. 



The research described in this paper was carried out at the Jet Propulsion Laboratory, 
California Institute of Technology, under a contract with the National Aeronautics and Space 
Administration. 
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Table 1. Determination of the centroid velocity of the Sun 



component 


n < 10 mas 


all stars 


Vx=e\ 


9.9 it 0.2 


10.5 ±0.1 


Vy = 


15.6 ±0.2 


18.5 ±0.1 


Vz = -el 


6.9 ±0.2 


7.3 ±0.1 
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Table 2. Determination of Ogorodnikov-Milne and higher degree parameters of the local 

velocity field 







n < 10 mas 


all stars 


Ll3 ■- 




5.91 ± 1.02 (5.8) 


6.21 ± 0.94 (6.6) 


L\2 '■ 


= B = h° 


-12.36 ± 1.26 (9.8) 


-13.36 ± 1.16 (11.5) 


L23 - 


= -hl 


0.13 ± 0.97 (0.1) 


-0.36 ± 0.89 (0.4) 


M12 


= A = 6e^^ 


14.08 ± 1.56 (9.2) 


13.83 ± 1.42 (9.8) 


M23 


= -3 62-1 


0.38 ± 1.37 (0.3) 


0.76 ± 1.25 (0.6) 


X2 = 


= X = -3e^ 


-0.63 ± 1.98 (0.3) 


1.02 ± 1.81 (0.6) 


Ml3 


= -3e^ 


-0.32 ± 1.38 (0.2) 


-2.13 ± 1.27 (1.7) 


= 


= C = 6e| 


-4.72 ± 1.64 (3.1) 


-3.03 ±1.43 (2.1) 


Other significant 


parameters 








-0.81 ± 1.28 (0.6) 


11.25 ±1.17 (9.6) 






-1.69 ±0.44 (3.8) 


-1.20 ±0.41 (2.9) 






0.52 ±0.24 (2.2) 


0.80 ±0.22 (3.7) 


el 




-0.075 ±0.041 (1.9) 


-0.100 ±0.037 (2.7) 



Note. — All parameters and their formal standard errors are spec- 
ified in kms" 1 kpc^^; the signal-to-noisc ratio is given in brackets. 
A solar velocity Vq = (9.9, 15.6, 6.9) kms-^ was subtracted for both 
sets of stars. 
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Fig. 1. — Dependence of mean orbital eccentricity on stellar age. Data extracted fr om the 
Geneva- Copenhagen spectroscopic survey of Hipparcos stars (INordstrom et al.ll2004l ). 
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Fig. 2. — The velocity field of Hipparcos stars generated by the vertical gradient of rotational 
velocity. 
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Fig. 3. — The velocity field of Hipparcos stars generated by the two unexpected magnetic 
vector harmonics, iff ^ and i?2"^- 



